B5.2 Applied PDEs 
Problem Sheet 2 
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Caustic - circle 
The boundary curve, solving for Po, Go 


x0 = Cos[s]; 


yO = Sin[s]; 
u0 = Cos[s]; 
S= 


Simplify[Solve[{D[u0, s] == pO0*D[x0, s] +q0*D[y0, s], p0*2+q0%*2 =: 1}, {p0, q0}]] 
{{p0 > 1, q0>0}, {p0 > -Cos[2s], q0 > -2Cos[s] Sin[s]}} 


Ray equations 


x=x0+2pO«t/. S[[2]]; 

y=y0+2q0x*t/. S[[2]]; 

Computing the envelope 

S2 = Simplify [Solve[Det[{{D[x, s], D[x, t]}, {D[y, s], D[y, t]}}] = 0, t]][[1]] 
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Simplify[x /. S82] 


1 

-—Cos[s] (-2+Cos[2s]) 
2 

Simplify[y /. 82] 


Sin[s]? 


Show[ParametricPlot[{x, y} /. S2, {s, 0, 27}, PlotStyle > Red], 
ParametricPlot[{x0, yO}, {s, 0, 27}, PlotStyle > Blue], PlotRange > All] 


